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Abstract 

It is shown that sufficiently smooth initial data for the Einstein- 
dust or the Einstein-Maxwell-dust equations with non-negative den- 
sity of compact support develop into solutions representing isolated 
bodies in the sense that the matter field has spatially compact sup- 
port and is embedded in an exterior vacuum solution. 

PACS: 04.20x 



1 Introduction 

While a detailed understanding of the motion of massive isolated bodies in 
general relativity is, for obvious physical reasons, of great interest, only a few 
special situations have been successfully analysed so far. 

The initial value problem for the Einstein equations with suitable matter 
fields, such as a perfect fluid model with an appropriate equation of state, 
is well understood in domains where the total energy density is positive. 
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Domains with vacuum-matter interfaces, which define the boundaries of the 
bodies, still pose technical problems. While quite general initial data for 
the Einstein-Euler system have been constructed by Dain and Nagy (0), 
which describe compact fluid balls embedded in an exterior vacuum field, 
the evolution problem for these data has not been solved yet in the same 
generality. 

If the density has a positive one-sided limit at the matter boundaries, 
where the pressure is required to vanish, the right hand side of Einstein's 
field equations will have discontinuities along the boundaries. Under the 
assumption of spherical symmetry, which is very special from the PDE point 
of view, the evolution problem has successfully been analysed in the presence 
of discontinuities by Ehlers and Kind ( |5j). Nothing is known for situations 
with lower degrees of symmetry. 

Situations where the density goes to zero at the boundaries with a certain 
smoothness have been discussed by Rendall (0). While no symmetries are 
assumed, his method requires the use of particular equations of state. These 
allow him to construct smooth solutions to the Einstein-Euler system with 
fluid bodies of spatially compact support, but, as pointed out by the author 
himself, the method he uses has undesirable properties such as a lack of 
uniqueness and the impossibility to discuss perturbations of static spherically 
symmetric solutions to the Einstein-Euler system. 

The present article deals with a further case in which the density goes 
to zero at the boundaries of the bodies and where no space-time symmetries 
are required. It assumes pressure free matter (also referred to as 'dust') or 
charged pressure free matter as a matter model. 

Two methods are known to write the Einstein-dust equations in hyper- 
bolic form. In 0] have been employed the frame formalism and the Bianchi 
equations for the conformal Weyl tensor to combine the hyperbolicity of the 
equations with a Lagrangian representation of the fluid flow. This is clearly 
of interest if one wants to control the location of the boundaries under gen- 
eral assumptions on the fluid model. These equations appear to be useful, 
however, only in four space-time dimensions. 

The representation of pQ starts from the Einstein equations in wave ('har- 
monic' in the older terminology) coordinates. To render the complete system 
hyperbolic, a further derivative is applied to the equations for the metric co- 
efficients so that one ends up again with a system of third order in the metric. 
The equations so obtained are applicable in all space-time dimension. Since 
the location of the boundary is not an essential problem in the situation 



2 



considered here these equations will be used in the following. 

In the existence proof we find it technically convenient to use a 'non- 
physical' extension of the matter flow vector field into the vacuum part of 
the solution. It is then shown that the solutions are in fact independent of 
these extensions and the latter can be ignored. Our main result concerning 
the motion of general relativistic bodies consisting of pressure free matter is 
then stated in the Theorem 12.31 The analogues result concerning charged 
pressure free matter is given in Theorem 13.31 

2 Solutions to the Einstein-dust equations 
2.1 The equations 

The stress energy tensor T of pressure-free matter is 

T a p = ruaup (2.1) 

where r is the scalar matter density, and u is the (future directed) flow vector 
field which satisfies 

u a u a = —1, hence u a SJ pu a = 0, (2.2) 

with V the covariant differential in the spacetime metric g. 
The tensor T satisfies the conservation law 

V a T aP = 0. (2.3) 

Its projection orthogonal to u, which reads ru a V a U[3 = 0, implies near points 
where r^O the geodesic equation 

u a V a u p = 0. (2.4) 

Its projection into the direction of u gives the continuity equation 

V a (ru a ) = u a d a r + rV a u a = 0. (2.5) 

The Einstein equations with pressure-free matter source read in n + 1 dimen- 
sional spacetime 

R a p = r{u a Ufs H T9ap)- (2-6) 

n — 1 

We shall regard equations (|2.4|) . ()2.5|) . ()2.6|) as our basic system of differential 
equations. 
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2.2 An existence and uniqueness theorem. 

The geometric initial data for the spacetime metric g on an initial manifold 
M are a Riemannian metric g and a symmetric 2-tensor k. The initial data 
for a dust source are a scalar function f on M and a tangent vector field 
v to M. A solution (V,g,r,u) of the coupled Einstein-dust equations is an 
Einsteinian development of the initial data set (M, g, k, f, v) if M can be 
diffeomorphically identified with an embedded submanifold of (V, g) so that 
g and k are respectively the induced metric and second fundamental form on 
M, while f is the function induced by r on M and v is the value on M of the 
dust velocity with respect to the proper frame (and the proper time) of an 
observer with timelike vector orthogonal 1 to M in (V, g). In local coordinates 
such that the values on M of the shift and the lapse of the development are 
respectively (3 = and N — 1, it holds that v l = (u°Y l u l , where u a are the 
components of u in the considered coordinate system at points of M. 

If M is compact, we denote by H s a usual Sobolev (Hilbert) space of func- 
tions or tensor fields on M which are square integrable together with their 
derivatives of order up to s in a given smooth Riemannian metric on the 
manifold M. The notation M s stands for continuous and bounded Rieman- 
nian metrics with derivatives in H s -\. On non compact manifolds suitable 
variants of these spaces can be used. Since the equations we consider are 
hyperbolic, the evolution problem can be localized. The discussion of the 
basic problem, namely the analysis of the evolution in domains containing 
boundaries of the bodies, will thus be essentially the same for compact or 
non-compact initial manifolds 2 . 

A tensor field of degree p on a manifold M x R, with R parametrized by 
t, can be decomposed, by projections on R or M, into a set ot t— dependent 
tensor fields, of degree 0, 1, . . . ,p on M. For a t-dependent a tensor field on 
M, t e [0,T[, the space E S (T) is defined as follows: 

E S (T) = C°([0, T[, H 9 ) n C\[0, T[, Hs-!) n C 2 ([0, T[, H s , 2 ). 

We say that a tensor on M x [0, T[ is in E S (T) if each of its projections is in 
E S (T). We denote by L S (T), respectively by U S (T), the space of Lorentzian 

1 Two such timelike vectors corresponding to two developments isometric under a dif- 
feomorphism f are mapped onto each other by f. 

2 Onc can replace, for instance, the spaces H s by spaces H" Joc of tensors with belong 
to H s in the open sets of a locally finite covering of the manifold M, with uniformly 
bounded H s norms. 
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metrics, respectively the space of unit vectors in the metric g, which belong 
to E S (T). 

We will prove the following local existence and geometrical global unique- 
ness theorem. 

Theorem 2.1 The Einstein equations in wave gauge with source a pressure 
free matter, form a hyperbolic Leray system for g, u and r, which is causal 
as long as g is Lorentzian and u is timelike. There is an interval [0, T[c R 
such that the Cauchy problem for these equations with data on the manifold 
M with g G M s , k G H s _i, v G H s _ x , f G H s _ 2 , where s > f + 2, \v\ g < I, 
has one and only one solution g G L S (T), u G Z7 a _i(T) ; r G E S _ 2 {T). 

Corollary 2.2 If the initial data satisfy the Einstein constraints, the solu- 
tion obtained in wave gauge satisfy the original Einstein-dust system. The 
solution is globally hyperbolic. There is a unique solution, up to isometry, in 
the class of maximal globally hyperbolic spacetimes if s > f + 3. 

Proof. The Einstein equations with dust source read in wave coordi- 
nates: 

R ap = -od^^a? + H aP {g, dg) = p afi = r(u a u p + ——g a p)- (2.7) 

They are to be coupled with equations (|2.4jl . (|2.5jl . 

We think of M as being embedded into the solution manifold and the 
coordinate x° to be chosen such that M = {x° = 0}. The solution manifold 
then takes close to M the form x M, with x° inducing a coordinate on 
the first factor, and M ~ {0} x M. The initial values for the unknowns 
g a p and u a are deduced from the geometric initial data set by assuming the 
contracted Christoffel symbols T a = g^" 1 Tp a 7 to vanish on M and the wave 
coordinates to be chosen such that the values on M of the lapse and the shift 
of the metric are respectively N — 1 and (3 = 0. 

We take a derivative of equations (|2.7j) in the direction of u and use 
equations (J2.4j) . (|2.5j) to obtain equations which are of third order in g but 
still do not contain derivatives of r: 

« 7 V 7 #S = ~r (u aUf3 + -l_^)V 7 u 7 . (2.8) 

The initial values for the second derivatives of g are determined on M so that 
equations (J2.7|) will hold on M. This will imply that the derivatives of the 
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contracted Christoffel symbols T a calculated from the derivatives of g will 
also vanish on M. 

To the equations (referred to below by their equation numbers) we assign 
the Leray-Volevic indices 

mfl22l = 0, mflHU) = 1, mO = 0, 

and to the unknowns g, u, r the indices 

£(g)=3, £(u) = 2, £(r) = 1. 

The matrix of principal parts of the various orders i — m is then diagonal 
and given by 

-\g^d^g aP \ 

u a d a u p 

u a d a r J 

The system is quasidiagonal with characteristic cotangent cone 

g^U^ta = o, 

the union of the light cone and of a spacelike hyperplane exterior to it, if u is a 
timelike vector. The system is thus hyperbolic and causal. The coefficients of 
the principal terms are in C 1 under the given hypothesis. The Leray-Dionne 
theory [3], [H] gives, after some work 3 , the existence of a number T > 
such that the considered Cauchy problem has a unique solution g e L S (T), 
ueU s ^(T),reE s _ 2 (T). 

Since equation (|2.8|) can then be rewritten as an ODE for w^l — p a p along 
the integral curves of u and we have (|2.7|) arranged to be satisfied on M, it 
follows that (J2.7j) is satisfied on the solution manifold. 

The proof of the corollary follows the same lines as in the vacuum case. 



2.3 The motion of isolated bodies. 

We have just proven that the Cauchy problem with data g, k, v, f as in The- 
orem EH is geometrically well posed and provides a unique solution g, u, r 
on V = [0,T[xM. There remains, however, an open problem. 

3 The lowering of the regularity of the data required by the general Leray theory of 
hyperbolic systems, which was worked out by Leray's student Dionne, was written up 
only in the case of one equation, though of arbitrary order. 
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Denote by uj the maximal open subset of M (assumed to be non-empty) 
on which f > 0. Suppose that its closure satisfies uj = supp(f) and, that f 
vanishes on an open set so that M \ uj 7^ 0. We can think then of w as a 
union of disjoint compact subsets of M which represent the space occupied by 
material bodies whose density tends continuously to zero at the boundary, 
since f G H 2 . While v represents the physically well defined 'matter flow 
velocity' at points of uj, it has no physical meaning at points where f vanishes 
(except, perhaps, on the boundary of uj as limit of physical velocities). In 
fact, the H s -i extension of v from uj to all of M has only be introduced as a 
convenient device to get the existence result. 

Equation (J2.5|) implies that r = along geodesies which start with tangent 
vectors u(p) with p G M \ uj. Along the geodesies with tangent vectors u{p) 
with p G uj the function r will be positive and bounded as long as the 
convergence -V a u Q remains bounded. If -V a M Q — ► oo, which indicates the 
development of a caustic of the geodesic vector field, we can expect r — ► oo 
and the evolution comes to an end. For geometric reasons we cannot have 
V a u a —>■ oo, whence r — > 0, in the future development. 

The closure of the 'geodesic tube' Q of points in V swept out by the 
geodesies with origin in uj coincides with the support of r in V. It represents 
the history of the material bodies. In V \ £1 the metric g satisfies Einstein's 
vacuum field equations. 

Theorem 2.3 Suppose f > 0. If s > ~ + 2 the physical solution in wave 
gauge obtained in Theorem \2.1\ which is given by the fields g and r onV and 
the field u on Q, is uniquely determined in a neighbourhood of M by the data 
g, k, f on M and the datum v on uj with uj = supp(f). It does not depend 
on the extension of v to M. 

Remark 2.4 The following discussion will also show that the life-time of the 
solution considered in Theorem VA. 1\ might be increased by suitable redefinitions 
of the field u outside Q. 

Proof. We note that the smoothness result of the theorem (ensured by 
choosing T small enough) excludes that any two of the geodesies correspond- 
ing to u cross each other on V (otherwise u would not even correspond to a 
well defined vector field). It follows in particular that none of the geodesies 
which start at points of M \ uj enters the set Q. Because Q is generated by 
time-like geodesies starting at points of uj, it has empty intersection with the 
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future g-domain of dependence D + (M \uj) in V of the set M \ uj. It follows 
from general results about the Einstein equations that g is determined on 
D + (M \ to) uniquely by the restriction of the data g and k to M \ uj. 

Suppose, is a datum on M analogues to v so that v* = v (whence 
«* = u) on uj and g,k,v*,r satisfy the requirements of Theorem 12.11 with 
s > | + 2. Then there exists a number T*, < < T, so that the tangent 
vectors of the ^-geodesies in V with tangent vector «* on M define a non- 
vanishing vector field u* on V* = [0,T,[xM C V. 

It follows that for points p G uj the restrictions to V* of the ^-geodesies 
with initial vectors it*(p) resp. tZ(p) coincide and are contained in V* fl Q. 
Equation (|2.5|) with u replaced by it* then determines a unique solution r* 
on V* which agrees with f on M. The functions r and r* have their support 
in the closure of V* fl Q, and we have in fact r* = r everywhere on V*. 

The expression on the right hand side of ()2.8|) remains unchanged if we 
replace u and r by and r*. The fields g, u*, r thus define on V* a solution 
to the Einstein-dust equations. By the uniqueness statement of Theorem 12. II 
it must coincide with the solution g*, it*, r* determined by the data g, k, iT*, f 
on M. Therefore g = g* and r t = r on 14, and, as seen previously u = it* on 
K fl IX The data g, k, v, f and g, k, v#, f determine the same physical solution 
on V*- 

The argument would be complete if we had T* = T for all possible exten- 
sions of v\ u to M. This can, however, not be expected. Since geodesic flows 
are defined by equations of second order, they tend to develop caustics. If 
this happens with the geodesies generating Q this might lead, as remarked 
above, to a blow up of the density r and to physical phenomena like shell 
crossings which have been studied in the literature (cf. 

Whether geodesies outside Q show this behaviour or whether geodesies 
starting on M \ uj tend to approach Q depends on the metric g as well as on 
the choice of extension of v\ w to M, and the latter will control to some extend 
the location of the set where this is going to happen. By a more judicious 
choice of the extension or by subsequent redefinitions of the extensions on 
suitable slices {t = const.} this phenomenon, which is of no physical or 
geometrical relevance but may restrict the domain of existence in a similar 
way as a bad choice of gauge, can be avoided. ■ 

The Einstein-dust constraints on M do not depend on the extension of v 
at points where f vanishes. They depend only of g, k, f on M and v on uj 
if supp(f) = uj. We can therefore state the following corollary. 
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Corollary 2.5 If the initial data g, k, f on M and v onto = supp(f) satisfy 
the Einstein - dust constraints, then the physically unique local solution ob- 
tained in wave gauge in theorem \2.1\ satisfies the full Einstein equations. It 
is locally geometrically unique. 

Proof. Since any extension to M of the given initial data satisfies the 
constraints, Corollary 12.21 shows that the corresponding solution on V ob- 
tained in wave gauge satisfies the full Einstein-dust equations. The same 
holds a fortiori for the restriction of this solution to the physically unique 
solution. 

To prove local geometric uniqueness we consider a solution of the Einstein- 
dust system g, r, u, in arbitrary coordinates, with Cauchy data g, k, f on M, 
and v on uj = supp(f) . One can show that there exist in a neighbourhood of 
M a change to harmonic coordinates reducing to the identity on M, hence 
preserving f and v, and such that g, k, are preserved. ■ 

3 Charged dust. 
3.1 The equations. 

The stress energy tensor of charged pure matter (dust) is the sum of the stress 
energy tensor of the matter and the Maxwell tensor of the electromagnetic 
field F : 

T a p = ru a up + T a /3, (3.1) 

with 

r a( 3 = F a x F px - ^g aP F^F^. (3.2) 
The Einstein equations read 

R a p = §af3 + r{u a u p + - g a p) (3.3) 

where we have set 

n — 3 1 

$ a/3 := F a x F/3\ + c n g aP F^F^, c n := (— - - - -). (3.4) 

The Maxwell equations are, with J the convection electric current of the 
density of charge q 

dF = and V.F = J, i.e. V Q F a/3 = J 13 := qu 13 ; (3.5) 
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they imply the conservation of charge equation 

\7 a (qu a ) = 0. (3.6) 

For simplicity we suppose the space-time to be simply connected so that 
there exists a 1-form A, the electromagnetic potential, with 

F = dA. (3.7) 

We take A in the Lorentz gauge, i.e. so that 

5 A = 0, i.e. V a A a = 0. (3.8) 

The Maxwell equations read then as a wave equation for A, namely 

V a V a A p - R\A X = J P = quP, (3.9) 

where we can replace the Ricci tensor R a p by the right hand side of ()3.3|) . 

Modulo the Maxwell equations (and u a u a = — 1), the stress energy con- 
servation equations are equivalent to 

V a (ru a ) = (3.10) 

and 

ru a V a vP + qu x F px = 0. (3.11) 
Near points where r ^ equations (J3.6)) and (|3.10|) imply 



u a d a (J) = 0, 



so that q/r is constant along the flow lines. It will be constant on the space- 
time to be constructed if it is constant initially. We will make this simplifying 
(though not necessary) assumption, and set 

q = cr, (3.12) 

with c some given number. Equation (|3. 1 lj) then implies near points where 
r > holds 

u^au 13 + cu x F /3X = 0. (3.13) 
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3.2 Existence and uniqueness theorem. 

Theorem 3.1 The Einstein equations in wave gauge with sources an electro- 
magnetic field of potential in Lorentz gauge, together with a charged pressure 
free matter, form a hyperbolic Leray system for g, A, q, u and r, which is 
causal as long as g is Lorentzian and u is timelike. There is an interval 
[0, T[c M such that the Cauchy problem for these equations with data on the 
manifold M such that g G M s , A G H s , k G H a -\, v G H s -i, f G H s -i, 
where s > ~ + 2, \v\g < 1 and q = cf with c = const., has one and only one 
solution g G L 3 (T), A G H S {T) u G U s -x{T), q = kr G E S ^ 2 (T). 

Corollary 3.2 If the initial data satisfy the Einstein and Maxwell constraints, 
the solution in wave and Lorentz gauges satisfies the original Einstein- Maxwell- 
dust system. It is globally hyperbolic. There is a unique, up to isometries, 
maximal, globally hyperbolic solution if s > f + 3. 

Proof. As in the previous section we differentiate the Einstein equations 
in wave gauge in the direction of u. We obtain, using the equation ()3.10|) (we 
could replace m 7 V 7 «a by its valued taken from 1)3.13)1 . but it is not necessary 
for our application) 

^V 7 ^ = f^, (3.14) 

with 

fop := m 7 V 7 $ Q/3 - r{(u a up + ^—^9ap) V 7 u 7 + u l V^(u a u^)}. (3.15) 

We also differentiate equations (J3.9)) in the direction of u to obtain 

M 7 V 7 V a V°A /3 - f x A x - RfaVyAx = cri-uNyV? + m 7 V 7 m /3 }. (3.16) 

We shall consider now the system for g, A, u, and r consisting of equations 
1)3.14)1 . 1)3.16)1 . ()3.13j) . and ()3.1(Jj) . To the equations we assign the Leray- 
Volevic indices 

m fTIl = 0. m flHIIBD = 0. mOSZEHD = 1, m fcUm = 0. (3.17) 
and to the unknowns g, A, u, r the indices 

£(g) = 3, £(A) = 3, £(u) = 2, £(r) = 1. (3.18) 
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We see that the matrices of principal parts is again diagonal, with the same 
kind of terms in the diagonal as in the previous section. The existence and 
uniqueness theorem in wave and Lorentz gauge follows from the Leray-Dionne 
theory. 

The proof of the existence of a solution of the full Einstein-Maxwell-dust 
system when the initial data satisfy the Einstein and Maxwell constraints 
follows the same lines as in the vacuum Einstein-Maxwell case. Geometric 
global uniqueness is also proved similarly. ■ 

3.3 The motion of isolated bodies. 

Again we denote by uj the maximal open subset of M on which f > and 
assume that its closure satisfies uj = supp(f). If f vanishes on an open set 
so that M \uj ^ 0, we denote by Q the subset of V = [0,T[xM which is 
generated by the integral curves of u which start at points of to. The support 
of r and q is then given by the closure Cl of Q in V. It follows that g and A 
satisfy the matter free Einstein-Maxwell equations in V \ fi. 

Theorem 3.3 Suppose that f > 0. The physical solution g, A, r, q on V 
and u on Q obtained in wave and Lorentz gauges in theorem 4-1 with initial 
data, on M, g , k, A, f , q = cf and an extension to M of the datum v on u 
with uj = supp(f) is independent, in a neighbourhood of M, of the extension 
ofv to M. 

Proof. The argument is the same as the proof of Theorem 12. HI with 
the geodesies curves considered there replaced by the integral curves of the 
vector field u defined by equation (j3.13|) . ■ 

Corollary 3.4 // the initial data g, k, A, f, q — cf on M, suppf C uj and 
v on uj, satisfy the Einstein - Maxwell - dust constraints, then the physically 
unique local solution obtained in wave and Lorentz gauge satisfies the full 
Einstein- Maxwell - dust equations. It is locally geometrically unique. 

Proof. The arguments are analogous to the proof of the corollary 12.51 
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